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(^Jf) 1 We develop a relativistic model to describe the bound states of positive 

■ energy and negative energy in finite nuclei at the same time. Instead of 

searching for the negative-energy solution of the nucleon's Dirac equation, 

CN 



we solve the Dirac equations for the nucleon and the anti-nucleon simultane- 
ously. The single-particle energies of negative-energy nucleons are obtained 



, through changing the sign of the single-particle energies of positive-energy 

CO 

' anti-nucleons. The contributions of the Dirac sea to the source terms of 

OS ' 

OS ■ the meson fields are evaluated by means of the derivative expansion up to 

the leading derivative order for the one- meson loop and one- nucleon loop. 
After refitting the parameters of the model to the properties of spherical 
nuclei, the results of positive-energy sector are similar to that calculated 



within the commonly used relativistic mean field theory under the no-sea 
approximation. However, the bound levels of negative-energy nucleons vary 
drastically when the vacuum contributions are taken into account. It implies 
that the negative-energy spectra deserve a sensitive probe to the effective 
interactions in addition to the positive-energy spectra. 

PACS number(s): 21.60.-n; 21.10.-k 
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I. INTRODUCTION 



The extension of the periodic system has a long history in chemistry and physics. 
Particularly at the end of the sixties and the beginning of the seventies, the time when 
heavy-ion collisions had its advent, the idea of super-heavy elements emerged. Besides 
the extension of the periodic system to the islands of super-heavy nuclei within the proton 
Z- and neutron iV-axis, there exist two fascinating direction, i.e., to extend Z and N into 
the negative sector, Z and N, and into the multistrangeness dimension. The first idea 
was put forward in detail in Ref. M where a collective anti-matter production mechanism 
was proposed. The main physical considerations are the following: In the relativistic 
treatment of nuclear phenomena, one uses the Dirac equation to describe nucleons in 
nuclei 



Here S and V are the scalar and time-like vector potentials acting on nucleons in the 
nuclear medium. In the simplest version of relativistic mean field theory the scalar and 
vector interaction are mediated by the sigma- and omega-meson exchange, respectively. 
The Dirac equation has two solutions, i.e., the positive-energy solution and the negative- 
energy solution. For the positive-energy solution, the nucleon central potential U cen ~ 
V + S while for the negative-energy solution, due to the "G-parity" [], the central anti- 
matter potential becomes U cen ~ —V + S. Relativistic mean-field (RMF) calculations 
in nuclear matter predict V ~ 300 MeV and S ~ —350 MeV ||. In these calculations 
the vacuum contributions have been neglected and the definite values of V and S are 
parameter-dependent. If these two values for V and S are inserted into the formulae 



1 The "G-parity" is an internal symmetry of the exchanged mesons in strong interactions, which 
connects the NN potential U(NN) = J2 m V™ to the NN potential U(NN) = £ m G(m)V m . G 
is the combination of the isospin symmetry and the C-conjugation rule and is usually defined 
as G = C exp(— infy). In particular, G(a) = G(p) = 1 and G{oS) = G(tt) = —1. 



[a • p + (5{M N + S) + V] * = EV. 
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of the central potentials, one immediately obtains U cen ~ —50 MeV and U cen ~ —650 
MeV. Therefore, the bound states of negative energy might be much deeper than the 
bound states of positive energy, although the exact potential depth of the bound states of 
negative energy has not been experimentally verified up to now. Furthermore, the vector 
potential V increases linearly with increasing of density. In violent relativistic heavy-ion 
collisions, the density may become very large, p ~ 5 — 10po- The potential of nucleons 
of negative energy (i.e., the anti- nucleons), may become larger than twice the nucleon's 
free mass with increasing density. Then, nucleons may be spontaneously emitted ||. But 
nucleons will also be emitted due to the dynamics, i.e., due to the time dependent change 
of their orbitals, i.e., due to the Fourier frequencies of the time dependent potentials and 
complex scatterings in the violent compression processes of a heavy-ion collision, and also 
due to temperature built-up in such violent encounters. These can create a great number 
of nucleon holes (i.e., anti-nucleons) in bound states. Then, "clusters" of holes (anti- 
nucleons) are distributed over shell-model-like orbitals. They are bound by the meson 
fields created by the compressed matter. They are, so to speak, prepared collectively 
for fusion into anti-matter clusters. Again, due to the Fourier frequencies of the violent 
nucleus-nucleus dynamics during and after the encounter, these anti-matter clusters can 
be kicked into the anti-cluster continuum, i.e., the negative-energy continuum, and may 
escape collectively out of the reaction zone. Another possible mechanism for anti-matter 
cluster escaping from the reaction zone is due to the space-time dependent event by event 
fluctuations of iV-body phase space of anti-nucleons which is not a smooth distribution 
function as in the case of 1-body calculations. Such collective creation processes of anti- 
matter clusters have a high potential to increase the probability for the production of light 
anti-matter clusters than the standard processes of particle scattering and coalescence. 

To realize the above theoretical conjecture for the production of anti-nuclei, - and 
analogously also for multi-A, multi-A nuclei, one should answer at least three questions in 
a rigorous way both from the theoretical side and the experimental side: Firstly, how deep 
are the bound states of negative energy (i.e., the bound states of anti-nucleons) in nuclei? 
Secondly, how does the potential of anti-nucleons increase with the increase of density 
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during the violent relativistic heavy-ion collisions and what is the dynamical production 
procedure for the anti-matter clusters? Thirdly, can we really have stable or meta-stable 
(heavy) anti-nuclei? In this work, we build a model to deal with the first question, which 
might be extended to answer the third question as well in the future. 

It is well known that relativistic mean-field theory had great success in describing 
the ground states of finite nuclei, i.e., the bound states of positive energy 0-0]. In this 
model, one usually considers only the positive-energy solution of the Dirac equation and 
neglects the contributions from the vacuum (the so-called no sea approximation) . With the 
parameters obtained as providing the best-fit to the properties of spherical nuclei, such as 
NL1 dJ and TM1 || parameter sets, the model describes well the ground-state properties 
of spherical nuclei as well as the deformation properties of even-even nuclei. Extensive 
investigation of various parameter sets in connection with predicting shells for superheavy 
nuclei have recently been given in Ref. 0. However, this all has been done in the no- 
sea approximation and without vacuum corrections. The later, the effects of quantum 
corrections, i.e., the vacuum contributions and their effects on the bound states of positive 
energy were investigated by several authors [|10|-|T4|| within the local density approximation 
for the one-nucleon loop, the one-meson loop and for the derivative expansion for the one- 
nucleon loop only (in a chiral model, loop corrections have been investigated in nuclear 
matter e.g., in Ref. [[HJ). The relativistic mean-field approximation is extended to the 
relativistic Hartree approximation (RHA). In that version, i.e., RHA, the parameters are 
fitted to the saturation properties of nuclear matter as well as the rms charge radius in 
40 Ca. The best-fit procedure within RHA to the properties of spherical nuclei has not been 
performed yet. In these preliminary studies it is found that the vacuum contributions do 
not improve the systematics of nuclei over RMF but the scalar density is decreased in 
the interior of nuclei. It, in turn, will cause the decrease of the scalar field S because 
the scalar field is coupled to the scalar density. Correspondingly, the vector field V will 
also decrease since the value V + S is controlled by the saturation properties. Therefore, 
even if the vacuum corrections may not cause significant difference to the bound states 
of positive energy after refitting the parameters, they should have strong influence on the 
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bound states of negative energy which are sensitive to the sum of scalar and vector field 
—V + S rather than the cancellation V + S (remember: V is positive, S is negative!). 

In the present paper we will develop a relativistic Hartree approach for the bound 
states of positive energy and negative energy in finite nuclei. The wave functions of 
the nucleons with positive energy are used to calculate the contributions of the valence 
nucleons. In principle, one should also use the wave functions of the anti-nucleons to 
evaluate the contributions of the bound states emerging from the Dirac sea. Considering 
the difficulties caused by the large number of anti-nucleons in bound states in any practical 
numerical procedure, we will employ the technique of derivative expansion developed in 



Ref. [16] to compute the effects of the Dirac sea, including the contributions of the one- 
nucleon loop as well as the one-meson loop. The paper is organized as follows: In Sect. 
II we introduce the effective Lagrangian used here. In Sect. Ill the plane-wave solutions 
of the Dirac equation in homogeneous nuclear matter are discussed. In Sect. IV we 
constitute the RHA approach for the bound states of positive energy and negative energy. 
In Sect. V we present the numerical results and discussions. Finally, a summary and 
outlook is given in Sect. VI. 

II. EFFECTIVE LAGRANGIAN 

We start from the Lagrangian density for nucleons interacting through the exchange 
of mesons 



£ = £ F + £i. (2) 

Here Cf is the Lagrangian density for free nucleon, mesons and photon 

£ F = #7^ - M N ty + Id^&o - U(a) - X -u)^ v 

+\mlu^ - + ±mjR„ ■ R" - \A, V A^ (3) 



and U(cr) is the self-interaction part of the scalar field [[17 



U(a) = \rn\a 2 + ^ba 3 + -ca*. (4) 
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In the above expressions ip is the Dirac spinor of the nucleon; a, u^, R M and A p represent 
the scalar meson, vector meson, isovector- vector meson field and the electromagnetic field, 
respectively. Here the field tensors for the omega, rho and photon are given in terms of 
their potentials by 

Unv = d^uj u - d u u^ (5) 
R-/ii/ = dnR v — c^R^, (6) 
A pu = d p A u — d u A p . (7) 

Ci is the interaction Lagrangian density 

1 - 1 - 

& = g^a - g^^ - -gp^i^r ■ ^ - -eip(l + Toh^A**. (8) 

Here r is the isospin operator of the nucleon and r is its third component. g CT , g w , g p 
and e 2 /47r = 1/137 are the coupling strengths for the a-, u-, p-meson and for the photon, 
respectively. M N is the free nucleon mass and m a , m W) m p are the masses of the a-, 
ou-, and p-meson. For simplicity, in the following discussions we consider the a- and u- 
exchange explicitly. The relevant formulae for the p-exchange and the electromagnetic 
field can be obtained in a straightforward way after calculating the tu-exchange. The 
detailed expressions for all these meson fields and the electromagnetic field will be given 
in Sect. IV. 

III. PLANE- WAVE SOLUTIONS OF THE DIRAC EQUATION IN STATIC 

NUCLEAR MATTER 

From the effective Lagrangian given in the above section, we obtain the Dirac equation 

[il^ - gu,!^ -M N + g a a} i> = 0. (9) 

In static nuclear matter, this reads 

d 

= H« • v + &M> + P(M N - g a a)} V>. (10) 
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We know that Eq. (10) has two solutions, i.e., the positive-energy solution E + and the 
negative-energy solution £L, which satisfy the following equations 



(£+-g^o)W(p) = [a-p + P(M N -g a a))U(p), 
(£L - g^ ) V(-p) = [a-p + /3{M N - gtr a)] V(-p). 



:n) 

12) 



Here U(p) and V(p) are two spinors |T3 



( 



U(p) = N 
where N is the normalization factor and 



\ ap v 



V(p) = iV 



£*(p)+m* 



(13) 



m* = M N - g a a, 



E*(p) = Jp 2 + m* 2 . 



Defining the effective positive-energy El and the effective negative-energy E*_ as 



(14) 
(15) 



E1 = E+- g^wo 



E*_ = E- - g^tuo, 



(16) 



together with Eqs. (11) ~ (13) one can perform calculations similar to the free Dirac 
equation and obtains 



El = E*(p) 



E*_ = -E*(p) 



(17) 



Therefore, 



E + = V P 2 + m * 2 + g^o, 
E_ = 



p 2 _|_ m *2 _|_ 



(18) 
(19) 



The wave packet ^ can be expanded according to the plane waves of positive energy and 
negative energy as 



d 3 p 



m 



1/2 



(2tt) 3 / 2 V£*(p) 



2 [&p,MP, s)e ip ' x -^+* + d+ s V(-p, s)e ip x - i ^ i ] . (20) 
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The coefficient b PyS are the probability amplitudes for waves with positive energy, whereas 
d+ s are those for negative energy. Both the nucleons with positive energy and negative 
energy transport forward in time. According to the Dirac's hole theory, each particle 
has its partner anti-particle which transports backward in time. One can write down the 
corresponding equations for the anti-nucleons which have the opposite signs for the energy 
and momentum term compared to the eigenequations of the nucleons 

- (E + + g^ ) V(p) = [-a • p + (M N - g a <j)\ V(p), (21) 

- [e. + g^ ) W(-p) = [-a • p + f3 (M N - g a a)\ U{-p). (22) 

Here E+, E_ are the positive energy and negative energy of the anti-nucleon and V(p), 
U(—p) are their corresponding eigenfunctions, respectively. Again, we define the effective 
positive-energy E+ and the effective negative-energy E*_ of the anti-nucleon as 

El = E + + g^ El = E_ + g^ . (23) 

With the spinors of Eq. (13) we obtain 

E* + = E*(p) El = -E*(p) (24) 

and 

E + = \Jp 2 + m* 2 - g^ujQ, (25) 

E_ = -^ p 2 + m *2 _ g ^ o . (26) 

From Eqs. (18), (19) and (25) , (26) one can easily find the relations 

E + = -E_ E_ = -E+. (27) 

Thus, in order to obtain the full spectrum of the Dirac equation, one can solve the 
equations for the nucleon with positive energy and negative energy, i.e., Eqs. (11) and 
(12). Alternatively, one can solve the Dirac equations for the nucleon and the anti-nucleon 
but both of them with positive energy, i.e., Eqs. (11) and (21). The wave packet ip can 
also be expanded according to the plane waves of the nucleon and the anti-nucleon as 
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#M)=/ (^(^yj ' EM(P.*)e*' x -^ t + <.V(p, a ) e - i ^ + ^*]. (28) 

Here now both the nucleoli as well as the anti-nucleon have positive energy. The nucleon 
transports forward in time while the anti-nucleon transports backward in time which can 
be seen from the sign of time in exponential functions. In relativistic quantum field theory, 
and d+ are explained as the creation operators for the nucleon and the anti-nucleon, 
respectively. This is the main strategy which will be used in our following considerations 
for finite nuclei. Instead of searching for two solutions with positive energy and negative 
energy for the Dirac equation of the nucleon, we solve the Dirac equations for the nucleon 
and the anti-nucleon simultaneously. For each equation we look for only one solution, 
that is, the positive-energy solution. 

When the negative energy of the nucleon is larger than the nucleon free mass, 
the system becomes unstable with respect to the nucleon-anti-nucleon pair creation. At 
zero-momentum, one has the critical condition 



1/2 



that is, 



if one defines 



g a a + guUo = 2M N , (29) 



V -S = 2M N (30) 



S = -g a cr V = g^ - (31) 



IV. RELATIVISTIC HARTREE APPROXIMATION OF FINITE NUCLEI 

In finite nuclei the meson fields in Eq. (10) are space-dependent. The field operator 
of the Dirac equation can be written as 

V(x, t)=E [WaWe"*' + tt(x)e ll » ( ] . (32) 

a 

Here the label a denotes the full set of single-particle quantum numbers. ^ a (x) are the 
wave functions of nucleons and ^(x) are those of anti-nucleons; E a and E a are their 
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positive energies, respectively. 6+ and d+ are nucleon and anti-nucleon creation operators 
that satisfy the standard anticommutation relations. We assume that the meson fields 
depend only on the radius and discuss the problem in spherically symmetric nuclei. In this 
case, the usual angular momentum and parity are good quantum numbers. As described 
in Refs. [|H||T!|], eigenf unctions of the angular momentum and the parity operator are the 
well-known spherical spinors. We make the following ansatz for the wave functions of 
nucleons 



(33) 



and anti-nucleons 



<(x) 



/ F a (r) (T t q 



jlm\ 



Here fL; m are the spherical spinors defined as WE 



(34) 



n jim = I m ' m sm) Y lm ,xi ma , (35) 

m'm s 

Yi m > are the spherical harmonics and Xi ma are t ne eigenfunctions of the spin operators. 
G a , F a and F a , G a are the remaining real radial wave functions of nucleons and anti- 
nucleons for upper and lower components, respectively. Applying the parity operator 
P = e l ^(3P , P changes x into — x, to Eqs. (33) and (34), one can easily find that the 
V'a(x) has the opposite eigenvalue of parity, (—1)', to the ^(x), (— 1)' +1 , as it should be 

II- 

Inserting Eqs. (33) and (34) into Eq. (10) (a — > a{r) and uj® — > u (r)), we obtain the 
coupled radial wave functions for the nucleon 

E a G a (r) = (--jL + ^ F a (r) + (M N - g ff a(r) + guJ uj (r)) G a (r), (36) 

E a F a (r) = (Jt + M G a (r) + (-M N + g(7 a(r) + guj u (r)) F a (r) (37) 

and the anti-nucleon 
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-E a F a (r) = (J- + G a (r) + (M N - g a <j{r) + &M>(r)) F a {r), (38) 
■4G a (r) = (-^ + F a (r) + + g.a(r) + ga M,(r)) G a (r) (39) 



respectively, where 



- ( i + i) fai = i +i _ (4Q) 

/ for j = / — | 



In order to resemble Schrodinger equations for the Dirac equations, we eliminate the small 
components. For the nucleon we eliminate the lower component while for the anti-nucleon 
the upper component. By defining the Schrodinger equivalent effective mass and potential 
of the nucleon 

M eff = E a -g w uj {r) + M N -g a (T{r), (41) 
U eff = M N - g a a(r) + g„cu (r) (42) 

and the anti-nucleon 

M eff = E a + gLJ LU (r) + M N -g a a(r), (43) 
U eff = M N - g a a(r) - g w uj (r), (44) 

we arrive at the Schrodinger equations for the upper component of the nucleon's wave 
function 



E a GJr) 



(~Tr + t) m <" {i + t) GJr) + u " lG " [r) (45) 



and the lower component of the anti-nucleon's wave function 

E a G a (r) = (-± + ^) M~ff (JL + ^) G a (r) + U eff G a (r). (46) 

The small components can be obtained through the following relations 

F a (r) = Mj f (± + !^)G a (r), (47) 
F a (r) = -M;/ f ^ + k -^G a (r). (48) 
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Of course, the radial wave functions are normalized 



poo 

/ dr 


1 G a (r) 


| 2 + 


Fa(r) \ 2 


Jo 








rod 








/ dr 


j G a (r) 


| 2 + 


F a (r) \ 2 


h 







(49) 
(50) 



From Eqs. (45) and (46) one finds that the Schrodinger equation of the anti-nucleon 
has the same form as that of the nucleon. The only difference relies on the definition of 
the effective mass and potential, that is, the vector field changes its sign. The so-called 
G-parity comes out automatically. Eqs. (45) and (46) can be solved numerically by the 
standard technique as described in Ref. ||. The single-particle energy of the nucleon and 
the anti-nucleon can be written as 



E a = r dr{G a (r) + M F a (r) + F a {r) (j- + M G a (r) 

Jo \ dr r J \dr r J 

+G a (r)U eff G a (r) - F a (r) {M eff - E a ) F a (r)}, 



dr{-G a (r) f -± + *2L ) F a (r) - F a (r) [± + ^ ] G a (r) 
o \ dr r J \dr r J 

+G a (r)U eff G a (r) - F a (r) (M e // - E a ) F a (r)}, 



(51) 



(52) 



which are obtained through the iteration procedure. The negative energies of nucleons 
are just the minus sign of E a . 

The meson fields in Eqs. (41 ) ~ (44) are determined by the Laplace equations 



V 2 - ml) a{r) = -g a p s (r) + ha 2 + ^ca 3 , 



ml) o> (r) 



;M r J 



(53) 
(54) 



where 



Ps{r) 



Po{r) 



Wi^iA + CT, 

=— oo 

Witf^ + CT, 



(55) 
(56) 



here Wi are the occupation numbers defined by the creation and annihilation operators 
of Eq. (32) as usual 0]. The CT are the counterterms. The sums in Eqs. (55) and 
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(56) run over the full spectrum of the Dirac equation: the negative-energy continuum 
(the positive-energy continuum of the anti-nucleon), some negative-energy bound states 
(positive-energy bound states of the anti-nucleon), the positive-energy bound states which 
correspond to the usual nucleon shell model states and the positive-energy continuum 
(included through the renormalization procedure). 

Let us first assume that the positive-energy continuum and the negative-energy con- 
tinuum together with the contributions of the counter terms will yield finite terms for 
Ap s (r) and Ap (r). The scalar and baryon density can then be expressed as 

Ps(r) = 7-^£w«(2j a + 1)(G» -Fl{r)) 

+^E^(% + 1)(<3J(0 - + Ap s (r), (57) 

Po(r) = + + F 2 a (r)) 



1 



X>/,(2j„ + l){G%r) + F$(r)) - Ap (r), (5* 



4irr 2 a 

The first terms on the right-hand-side of Eqs. (57) and (58) denote the contributions of 
the shell model states while the second terms represent the contributions of the bound 
states of negative energy. If one thinks about that there might exist twenty thousand 
nucleons in the bound states of negative energy, it is obviously not practical to compute the 
contributions of those bound states through evaluating their wave functions. Fortunately, 



a rather elegant technique has been developed by several authors [16|], which takes into 
account the vacuum contributions to the source term of the meson fields in finite nuclei. 
Let us write 

Ps(r)=p v s al (r)+p s s ea (r), (59) 
Po(r)=p v al (r)+p s ea (r), (60) 

where p v s al (r) and pQ al (r) are just the first terms on the right-hand-side of Eqs. (57) and 
(58). p s g a (r) and p S Q a (r) are the contributions of the vacuum, including the bound states 
and the continuum as well as the counterterm contributions. Originally, vacuum correc- 
tions for finite nuclei (in one-loop approximation) were included only in a local density 
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approximation . Later on, Perry Jy]] and Wasson [|T^] considered derivative corrections 



to the nuclear matter results. It was found that the leading derivative correction is of the 
same order of magnitude as the effective potential (i.e., the nuclear matter results) while 
the next-order derivative correction is two or three orders of magnitude smaller than the 
leading order. That means that the derivative expansion converges rapidly. Since we 
include the non-linear self-interaction of the scalar field in the model, the contributions 
of the one-meson loop from the scalar meson should be taken into account in addition to 
the one-nucleon loop. The contributions of the one-meson loop have been calculated in 



Ref. JT3| up to the effective potential term. Here this is extended to include the leading 
derivative correction. 

The effective action of the system at the one-loop level can be written as 

T = Jd A x (^dpaVa - U(a) - ~cu^ + -m*^" + CT) + r valcncc + T«(a) + r<%). 

(61) 

Here r va i ence is the contribution from the valence nucleons, which for time independent 
background fields is just minus the energy of the valence nucleons. T^(a) and [jp) 
represent the contributions of the Dirac sea stemming from the one-meson loop and one- 
nucleon loop, respectively. They are defined as |?U| 

i 



r«((7) = -ftTrln 
v ; 2 



P 2 - (ml + ba + ^ca 2 ) 



(62) 

r« = -ihTr In [P - m* - g u </] , (63) 



where the trace is over spatial and internal variables. r^'(a) and r^ 1 -* (ip) can be expanded 
in powers of the derivatives of the scalar and vector fields. Using Lorentz invariance one 
can determine the functional Taylor series as 

T ( V) = / d'x (-if V) + VV)(<V) 2 + r«(a)(<Vx) 4 + , (64) 
r«(V0 = / d'x (-Vj'V) + ^V)0V) 2 + \z^(a)u^ + ...) . (65) 

Here v£\a) and (a) are the effective potentials from the one-boson loop and one- 
fermion loop, in which the field is a constant, u(x) = a , the same situation as in nuclear 
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matter. These two terms contain divergent part and should be renormalized. Through 
adding the suitable counterterms, vlp (o) and V F (o) can be evaluated in nuclear matter 



which turn out to be fl21 



m„ 



■57TJ' 

§( 

1 



bo co 2 
— r + 



bo 



m* a 2m 2 
2 \2 n 



' /\ fecr co- 
in 1 + — + 



ccr 



6a 



mi 



2ml 



mi 



mi 



bo 



mi 



2m\ 

3co 2 \ 
2ml I 



bo co 
— r + 



mi 



(M 



N 



4tt 2 
+ ^-M N gio 3 



l<jo) In ( 1 
25 



M: 



N 



+ M%g a o 



1 (bo 
12 [ril 



;M^g> 2 



2m 2 

41 



(66) 
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4 4 



(67) 



As discussed above, the derivative expansion converges rapidly. Thus, here we consider 
only the lowest order derivative terms in Eqs. (64) and (65). By means of the tech- 
nique developed in Ref. [16], one can determine the functional coefficients of the leading 
derivative correction, which read 

1 {b + co) 2 



1216vr 2 (m2 + bo + \co 2 )' 



2tt 2 

r 2 



2 

^ln 



77T 



M 



N 



In 



/ m* 



3tt 2 \M, 



N 



(68) 
(69) 
(70) 



The above expressions are finite, and independent of any renormalization conditions. The 
divergent integral appearing in T^(o) and T^(ip) are included in the effective potential 
terms of vjp(a) and V F l \o). The equations of meson fields, i.e., Eqs. (53) and (54), 
can be obtained through minimizing the effective action of Eq. (61) with respect to the 
corresponding fields. With the definitions of Eqs. (57) ~ (60), r va i encc contributes to 
the p v s al (r) and pQ al (r) while T^(o) and (ip) contribute to the p| ea (r) and p S Q a {r). 
Inserting Eqs. (64) ~ (70) into Eq. (61) and minimizing the effective action, one obtains 
the concrete expressions of p| ea (?") and p S Q a {r) which read 



]_d_ 

ga 00 



v£\o) + V F '\o 



(1), 



ldZ^(o) 
2 do 



yo) 2 - v • zM(o)Vo 
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and 

dZW(a) 1 
da ~ 192vr 2 

Note that Po ea (r) * s a total derivative and thus the baryon number is conserved. 

To include the contributions of the p-exchange and the electromagnetic field, in Eqs. 
(41) and (42) we make a replacement of 

gwW (r) — ► &M)(r) + ^g p T 0a R 0fl (r) + ^e(l + T 0a )A (r) 

and in Eqs. (43) and (44) 

gu>v (r) — > gu,^o(r) - ^g p T 0a R 0fi (r) + ^e(l + r 0a )A (r). 

Here we have defined that the anti-particle has the same isospin factor as the correspond- 
ing particle. Note that the G-parity of p-meson is positive. The field equations of the 
p-meson and the photon read 

(V 2 - ml)R 0fi (r) = -\g p {p v Ar) + p^(r)), (76) 
V%(r) = -e« (r) + p^ (r)). (77) 

We assume isospin-symmetry in the bound states of negative energy, therefore, p s qq{t) = 0. 
Other densities appearing in Eqs. (76) and (77) can be calculated with the same steps as 
given above. At the end we obtain 
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2c(b 



co) 



[mi 



bo 



[mi 



(b + c o) 
+ bo 



(75) 




(78) 



(79) 



(80) 



If one includes the p-exchange and the electromagnetic field in Eq. (63), a term 



should be added to Eq. (71), which represents the contribution of the electromagnetic field 
to the derivative correction of the one-nucleon loop. The corresponding contribution of the 
p-meson field is negligible in the isospin-symmetry vacuum. The Schrodinger equations 
(45), (46), and the Laplace equations (53), (54), (76), (77), must be solved numerically 
in a self-consistent iteration procedure [§]]. 



Since we employ the derivative expansion to evaluate the contributions of the Dirac 
sea to the source terms of the meson fields, the wave functions of anti-nucleons, which are 
used to calculate the single-particle energies, are not involved in evaluating the vacuum 
contributions to the scalar and baryon density which are, in turn, expressed by means of 
the scalar and vector field as well as their derivative terms. The Dirac equation of the 
nucleon and the equations of motion of mesons (containing the densities contributed from 
the vacuum) are solved within a self-consistent iteration procedure M. Then, the Dirac 
equation of the anti-nucleon is solved with the known mean fields to obtain the wave 
functions and the single-particle energies of anti-nucleons. The space of anti-nucleons 
are truncated by the specified principal and angular quantum numbers n and j with the 
guarantee that the calculated single-particle energies of anti-nucleons are converged when 
the truncated space is extended. We find that the results are insensitive to the exact 
values of n and j provided large enough numbers are given. We have used n = 4, j = 9 




V. NUMERICAL RESULTS AND DISCUSSIONS 



for 16 0; n = 5, j = 11 for 40 Ca; and n = 9, j = 19 for 208 Pb. 



17 



As pointed out in the Introduction, in the previous RHA calculations for the bound 
states of positive energy [Hill], the parameters of the model are fitted to the saturation 
properties of nuclear matter as well as the rms charge radius in 40 Ca. The best- fit routine 
within the RHA to the properties of spherical nuclei has not been performed yet. Thus, 
we first fit the parameters of the effective Lagrangian presented in Sect. II within the 
RHA to the empirical data of binding energy, surface thickness and diffraction radius of 
eight spherical nuclei 16 0, 40 Ca, 48 Ca, 58 Ni, 90 Zr, 116 Sn, 124 Sn, and 208 Pb as has been done 
in Ref. for the RMF model. The experimental values for the observables used in the 
fit are given in Table I. In the fitting processes, we distinguish two different cases with 
(RHAl) and without (RHAO) nonlinear self-interaction of the scalar field. The obtained 
parameters and the corresponding saturation properties are given in Table II. For the sake 
of comparison, two sets of the linear (LIN) and nonlinear (NL1) RMF parameters from 
Ref. 0] are also presented. One can see that the RHA gives a larger effective nucleon mass 
than the RMF does, which is mainly caused by the feedback of the vacuum to the meson 
fields, as can be seen from Eqs. (71) ~ (74). When the effective nucleon mass decreases, 
the scalar density originated from the Dirac sea p| ea increases. It is negative and cancels 
part of the scalar density contributed from the valence nucleons p v s al , which causes the 
effective nucleon mass to increase again. At the end, it reaches a balance value. In the 
fitting procedure, we have tried different initial values giving smaller effective nucleon 
mass. After running the code many times, all of them slowly converge to a large m* . 
The larger effective nucleon mass explains why a larger \ 2 value is obtained for the 
RHAl compared to the NL1. If one uses the current nonlinear RMF/RHA models to fit 
the ground states properties of spherical nuclei, an effective nucleon mass around 0.6 is 
preferred. The situation, however, might be changed when other physical ingredients, e.g., 
tensor coupling for the vector fields, correlation effects, three-body forces, are taken into 
account, which warrants further investigation. On the other hand, in the case of linear 
model, the RHAO gives a better fit than the LIN does. This is mainly due to the vacuum 
contributions which improve the theoretical results of the surface thickness substantially, 
and finally improve the total x 2 value. 
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Fig. 1 displays the equations of state of nuclear matter as well as the scalar and vector 
potentials as a function of density calculated with four sets of parameters given in Table 
II. It is clear that the RHA exhibits a softer equation of state compared to the RMF, 
mainly due to the larger m*. In the mean time, the strengths of the scalar and vector 
potentials decrease substantially in the RHA as expected. The same situation happens 
at finite nuclei. In Fig. 2 and Fig. 3 we present the scalar and vector potentials in 16 0, 
40 Ca and 208 Pb computed with the NL1 and the RHA1 set of parameters, respectively. 
The potentials calculated with the RHAl are about half of that calculated with the NL1, 
implying a strong feedback of the vacuum to the meson fields. 

As given in Sect. IV, the contribution of the Dirac sea to the baryon density is a total 
derivative. The net baryon number is conserved. Fig. 4 depicts the fractions of the baryon 
density stemming from the vacuum in spherical nuclei of 16 0, 40 Ca and 208 Pb. It can be 
seen that the p S Q a is more pronounced in light nuclei. In heavy nuclei it is negligible. 
Specifically, the p s ea /p < 4.0% in 16 0, 2.3% in 40 Ca and 0.6% in 208 Pb. In Fig. 5(a) 
we compare the baryon densities of 16 calculated within the RHA and the RMF model, 
respectively. One can see that the difference is not very siginificant. In the case of the 
RHA, the contributions of different sources to the baryon density are shown in Fig. 5(b). 
The vacuum contribution changes its sign from the interior to the surface of the nucleus. 
At large r, the Dirac-sea effect is negligible. This can be observed in Fig. 4 too. For 
different nuclei, after the r exceeds the typical values of surface range, the p S Q a decreases 
rapidly. 

In Fig. 6 the charge densities of three spherical nuclei computed with the NL1 and 
the RHAl set of parameters are compared with the experimental data. It seems that for 
light nucleus 16 the results of the NL1 are closer to the data than that of the RHAl. 
Alternatively, for media and heavy nuclei, the charge densities calculated with the RHAl 
show a better agreement with the data and less shell fluctuation in the interior of nuclei 
than that with the NL1 model. The shell fluctuation can be best expressed via the charge 
density in 208 Pb as 

5p = p c (1.8 fm)-p c (0.0 fm). (81) 
19 



The empirical value is —0.0023 fm -3 [|J], which is nicely reproduced in the RHA (see Table 
II) while the RMF overestimates Sp by a factor of 3, sharing the same disease with the 
non-relativistic mean field theory |22| . 

The effects of the Dirac sea on the physical quantity, specifically, the binding energy 
per nucleon are displayed in Table III for eight spherical nuclei. The calculations are 
performed with the RHA1 set of parameters. It can be found that the theoretical pre- 
dictions including the vacuum contributions are in good agreement with the empirical 
values. The magnitudes of the Dirac-sea corrections are quite similar for different nuclei. 
The absolute ratios of the Dirac-sea effects and the binding energy stay between 16% 
and 17%. All corrections decrease the binding energy per nucleon. The contributions of 
the derivative terms are in the same order of magnitude as the effective potential terms. 
When the nucleus becomes heavier the derivative terms turn out to be smaller since the 
fields become more stable. In all cases, the effects of the one-nucleon-loop is about four 
to five times of that of the one-meson-loop. 

Now let us go to the single-particle levels. In Table IV and V we present the re- 
sults of both positive- and negative-energy proton and neutron spectra of 16 0, 40 Ca and 
208 Pb. The binding energy per nucleon and the rms charge radius are given too. The 
numerical calculations are performed within two frameworks, i.e., the RHA including the 
contributions of the negative-energy sector to the source terms of the meson fields and 
the RMF taking into account only the valence nucleons as the meson-field sources. The 
experimental data are taken from Ref. ||23|| . From the table one can see that all four sets 
of parameters can reproduce the empirical values of the binding energies, the rms charge 
radii and the single-particle energies of positive-energy states fairly well. For the E/A and 
the r c /i, the agreement between the theoretical predictions and the experimental data are 
improved from the LIN to the RHAO, RHAl and NL1 set of parameters. For the spectra 
of positive-energy states, due to large error bars, it seems to be difficult to queue up the 
different sets of parameters. However, because of the large effective nucleon mass, in the 
current models the RHA has smaller spin-orbit splitting (see lpi/2 and IP3/2 state) com- 
pared to the RMF. This situation can be improved through introducing a tensor coupling 
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for the u> meson J| . With a suitable chosen coupling strength for the tensor term, a rea- 
sonable spin-orbit splitting may be obtained while a large m* remains. The effects of the 
tensor-coupling terms will be investigated in the future studies. For the negative-energy 
sector, no experimental data are available. In all four cases, the potentials of negative- 
energy nucleons are much deeper than the potentials of positive-energy nucleons. On the 
other hand, one can notice the drastic difference between the RHA and the RMF calcu- 
lations - the single-particle energies calculated from the RHA are about half of that from 
the RMF as can be expected from Fig. 2 and 3, exhibiting the importance of taking into 
account the Dirac sea effects. It demonstrates that the negative-energy spectra deserve 
a sensitive probe to the effective interactions in addition to the positive-energy spectra. 
The spin-orbit splitting of negative-energy sector is so small that one nearly can not dis- 
tinguish the lpi/2 and the 1^3/2 state. This is because the spin-orbit potential is related 
to d(S + V)/dr in the negative-energy sector and two fields cancel each other to a large 
extent. Nevertheless, the space between the Is and the lp state is still evident, especially 
for lighter nuclei. This might be helpful to separate the process of knocking out a lsi/2 
negative-energy nucleon from the background - a promising way to measure the potential 
of the anti-nucleon in laboratory. 

Fig. 7 depicts the potentials of the anti-nucleon in 208 Pb computed within the NL1 
and the RHAl model. The corresponding potentials of the nucleon are presented too. 
One can easily find that the NL1 and the RHAl set of parameters give the similar poten- 
tials for the nucleon (except in the center of the nucleus where the results of the RHAl 
model exhibit certain fluctuations caused by the p s g a contributed from the vacuum) while 
the potentials of the anti-nucleon differ substantially. A much weaker anti-nucleon po- 
tential is obtained in the RHAl model. This results in a large critical density (around 
9.5po) for the spontaneous nucleon-anti-nucleon pair creation as can be seen in Fig. 8 
where the density dependence of the nucleon and the anti-nucleon energy in symmetric 
nuclear matter is given. It should be mentioned that the results of Fig. 7 and 8 are very 
sensitive to the effective nucleon mass which can not be determined unambiguously in a 
model-independent way via experiments. In order to know the individual scalar and vec- 
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tor potentials, one has to analyse the empirical data both from the positive-energy and 
the negative-energy sector. The later is, unfortunately, currently unavailable. Further 
investigation is apparently needed before coming to a definite conclusion. 

VI. SUMMARY AND OUTLOOK 

This paper develops a model to describe the bound states of positive energy and 
negative energy in nuclei. Instead of expanding the field operator of the Dirac equation 
in plane waves of nucleons with positive energy and negative energy, here the expansion 
is done with plane waves for nucleons and anti-nucleons separately, both of them with 
positive energy. In this case, the Dirac equations of both, the nucleon and the anti-nucleon, 
can be reduced to Schrodinger-equivalent equations. For each equation, we then search for 
only one solution, i.e., the positive-energy solution. The numerical procedure is similar to 
the one currently used in the relativistic mean field theory for the bound states of positive 
energy, except that one more equation for the anti-nucleon is implemented. Thus, the 
model is solvable with existing techniques. The single-particle energy of the nucleon with 
negative energy is just the negative of the single-particle energy of the anti-nucleon with 
positive energy, i.e., the solution of the anti-nucleon's Schrodinger equation. 

The contributions of the Dirac sea to the source terms of meson fields can be sepa- 
rated into two parts, that is, the contributions of the negative-energy bound states and 
the negative-energy continuum. In principle, one should use the wave functions of anti- 
nucleons to evaluate the effects of the bound states emerging from the Dirac sea. There 
might exist a large number of anti-nucleons in bound states. Hence, it is apparently not 
practical to calculate all those wave functions in the iteration procedure. Here we employ 
the technique of derivative expansion for the one-meson and one-nucleon loop to take into 
account the vacuum contributions of both the bound states and the continuum. 

In the numerical procedure, we first fit the parameters of the model to the properties of 
spherical nuclei. Two sets of parameters with and without nonlinear self-interaction of the 
scalar field are obtained for the present RHA model including the vacuum contributions. 
They are then used to investigate the vacuum polarization effects on both positive-energy 
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and negative-energy sector. The corresponding calculations of the RMF model are also 
presented for comparison. Our results show that both the RMF and the RHA model 
describe the properties of spherical nuclei very well. Due to the feedback of the vacuum 
to the meson fields, the scalar and vector potentials decrease in the RHA. This causes 
the drastic difference on the single-particle energies of negative-energy nucleons calculated 
within the RHA model and within the RMF model, while the single-particle energies of 
positive-energy nucleons coincide each other within two models. Since the negative-energy 
sector is sensitive to the sum of the scalar and vector field — V + S while the positive- 
energy sector is sensitive to the cancellation of the fields V + S, the study of both of 
them in an unified framework will lead to the determination of the individual S and V\ 
Thus, it is currently very important to have experimental data to check the theoretical 
predicted bound levels of negative energy. It will provide us with a chance to judge the 
physical necessity of introducing strong scalar and vector potentials in the Dirac picture. 
If this picture is valid for the nucleon-nucleus and anti-nucleon-nucleus interactions, a 
fascinating direction of future studies is to investigate the vacuum correlation and the 
collective production of the anti-nuclei in relativistic heavy-ion collisions. Experimental 
efforts in this direction are presently underway P3 . 
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TABLE I. The experimental values for the observables included in the fit, the binding energy 
Eb, diffraction radius R and surface thickness a. In the last line we also give the adopted errors 
AO„ for the fit. 





E B (MeV) 


R (fm) 


a (fm) 




±00 


— 127.6 


2. .777 


0.839 




40 Ca 


-342.1 


3.845 


0.978 




48 Ca 


-416.0 


3.964 


0.881 




58 Ni 


-506.5 


4.356 


0.911 




90 Zr 


-783.9 


5.040 


0.957 




116 Sn 


-988.7 


5.537 


0.947 




124 Sn 


-1050.0 


5.640 


0.908 




208p b 


-1636.4 


6.806 


0.900 




AO n /O n 


0.2% 


0.5% 


1.5% 




TABLE II. Parameters of the RMF and the RHA models as well as the correspondin 


g satura- 


tion properties. The results of shell fluctuation 


and the x 2 


values of different sets of parameters 


are also presented. 












RMF 




RHA 






LIN 


NL1 


RHA0 


RHA1 


M N (MeV) 


938.000 


938.000 


938.000 


938.000 


m a (MeV) 


615.000 


492.250 


615.000 


458.000 


m w (MeV) 


1008.00 


795.359 


916.502 


816.508 


m p (MeV) 


763.000 


763.000 


763.000 


763.000 


ga 


12.3342 


10.1377 


9.9362 


7.1031 


goj 


17.6188 


13.2846 


11.8188 


8.8496 


gp 


10.3782 


9.9514 


10.0254 


10.2070 


b (fm" 1 ) 


0.0 


24.3448 


0.0 


24.0870 


c 


0.0 


-217.5876 


0.0 


-15.9936 


Po (fm~ 3 ) 


0.1525 


0.1518 


0.1513 


0.1524 


E/A (MeV) 


-17.03 


-16.43 


-17.39 


-16.98 


m*/M N 


0.533 


0.572 


0.725 


0.788 


K (MeV) 


580 


212 


480 


294 


a 4 (MeV) 


46.8 


43.6 


40.4 


40.4 


Sp in 208 Pb (fm" 3 ) 


-0.0075 


-0.0070 


-0.0016 


-0.0030 


x 2 


1773 


66 


1040 


812 
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TABLE III. Experimental and theoretical binding energy per nucleon as well as the vacuum 
corrections. The calculations are performed with the RHA1 set of parameters. The Dirac sea 
effects are separated by the effective potential terms and the derivative terms as well as the 
contributions of the one-nucleon-loop and the one-meson-loop. 



Exp. Theory Dirac Sea Eff. Pot. Deri. Terms Nucl. Loop Meson Loop 



16Q 


-7.98 


-8.00 


1.37 


0.58 


0.80 


1.10 


0.27 


40 Ca 


-8.55 


-8.73 


1.43 


0.82 


0.60 


1.17 


0.26 


48 Ca 


-8.67 


-8.51 


1.39 


0.82 


0.57 


1.13 


0.25 


58 Ni 


-8.73 


-8.44 


1.44 


0.87 


0.57 


1.18 


0.26 


90 Zr 


-8.71 


-8.74 


1.42 


0.96 


0.46 


1.16 


0.25 


116 Sn 


-8.52 


-8.61 


1.39 


0.98 


0.42 


1.14 


0.25 


124 Sn 


-8.47 


-8.50 


1.34 


0.96 


0.39 


1.10 


0.24 


208p b 


-7.87 


-7.93 


1.30 


0.98 


0.32 


1.07 


0.23 
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TABLE IV. The single-particle energies of both positive- and negative-energy protons as well 
as the binding energy per nucleon and the rms charge radius in 16 0, 40 Ca and 208 Pb. 





RMF 




RHA 








LIN 


NL1 


RHA0 


RHA1 


EXP. 


16Q 

E/A (MeV) 


7.80 


8.00 


8.01 


8.00 


7.98 


r c h (fm) 


2.59 


2.73 


2.62 


2.66 


2.74 


POS. ENE. 












ls 1/2 (MeV) 


42.99 


36.18 


32.21 


30.68 


40±8 


lp 3 /2 (MeV) 


20.71 


17.31 


16.09 


15.23 


18.4 


lp 1/2 (MeV) 


10.85 


11.32 


12.98 


13.24 


12.1 


NEG. ENE. 












ls V 2 (MeV) 


821.30 


674.11 


413.62 


299.42 




lp 3 /2 (MeV) 


754.62 


604.70 


369.78 


258.40 




lp 1/2 (MeV) 


755.43 


605.77 


370.36 


258.93 




40 Ca 












£7/4 (MeV) 


8.38 


8.58 


8.65 


8.73 


8.55 


reft (fm) 


3.36 


3.48 


3.39 


3.42 


3.45 


POS. ENE. 












ls l/2 (MeV) 


51.21 


46.86 


38.64 


36.58 


50±11 


lp 3/2 (MeV) 


35.05 


30.15 


27.11 


25.32 




l Pl/2 (MeV) 


29.25 


25.11 


25.17 


24.03 


34±6 


NEG. ENE. 












ls V 2 (MeV) 


840.76 


796.09 


456.58 


339.83 




lp 3 /2 (MeV) 


792.36 


706.36 


424.85 


309.24 




lp 1/2 (MeV) 


792.75 


707.86 


425.14 


309.52 




208p] 3 












£/4 (MeV) 


7.83 


7.89 


7.96 


7.93 


7.87 


(fm) 


5.34 


5.52 


5.43 


5.49 


5.50 


POS. ENE. 












lsi/2 (MeV) 


58.71 


50.41 


44.43 


40.80 




lp 3/2 (MeV) 


52.74 


44.45 


39.87 


36.45 




lp 1/2 (MeV) 


51.83 


43.75 


39.49 


36.21 




NEG. ENE. 












lsi/2 (MeV) 


830.16 


717.01 


476.61 


354.18 




lp 3 /2 (MeV) 


819.15 


705.20 


466.08 


344.48 




lp 1/2 (MeV) 


819.22 


705.28 


466.13 


344.52 
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TABLE V. The single-particle energies of both positive- and negative-energy neutrons. 





RMF 

LIN 


NL1 


RHA 

RHA0 


RHA1 


EXP. 


16Q 

POS. ENE. 

lsi /2 (MeV) 
lp 3/2 (MeV) 
l Pl/2 (MeV) 
NEG. ENE. 
lsi/2 (MeV) 
lp 3 /2 (MeV) 
lp 1/2 (MeV) 


47.23 
24.70 
14.74 

814.98 
748.48 
749.22 


40.21 
21.07 
15.01 

667.93 
598.74 
599.74 


36.33 
19.99 
16.86 

407.44 
363.83 
364.37 


34.71 
19.04 
17.05 

293.23 
252.48 
252.97 


45.7 
21.8 
15.7 


40 Ca 
POS. ENE. 

lsi/2 (MeV) 
lp 3 /2 (MeV) 
l Pl/2 (MeV) 
NEG. ENE. 
ls 1/2 (MeV) 
lp 3 /2 (MeV) 
lp 1/2 (MeV) 


59.36 
42.94 
37.17 

828.82 
781.18 
781.46 


54.85 
37.79 
32.73 

783.87 
694.80 
696.18 


46.67 
34.91 
32.99 

444.85 
413.71 
413.93 


44.48 
32.98 
31.71 

327.96 
298.04 
298.26 




208 pk 

POS. ENE. 

lsi/2 (MeV) 
lp 3 /2 (MeV) 
lp 1/2 (MeV) 
NEG. ENE. 
ls 1/2 (MeV) 
lp 3/2 (MeV) 
lp 1/2 (MeV) 


65.19 
58.50 
57.73 

789.37 
779.44 
779.45 


58.97 
52.44 
51.82 

678.23 
667.70 
667.73 


50.99 
46.05 
45.71 

435.30 
425.81 
425.82 


47.40 
42.66 
42.45 

313.18 
304.61 
304.61 
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FIGURE CAPTIONS 



Fig.l Equations of state of nuclear matter as well as the scalar and vector potentials 
calculated within the RMF and the RHA models under different sets of parameters 
given in Table II. 

Fig. 2 The scalar potentials in 16 0, 40 Ca and 208 Pb. 

Fig. 3 The vector potentials in 16 0, 40 Ca and 208 Pb. 

Fig. 4 The fractions of the baryon density contributed from the Dirac sea. The calcula- 
tions are performed for 16 0, 40 Ca and 208 Pb with the RHAl set of parameters. 

Fig. 5 The upper panel displays the baryon density in 16 computed with the NL1 and 
the RHAl set of parameters, respectively. The lower panel shows the contributions 
of the valence nucleons and the Dirac sea to the baryon density. The calculations 
are performed for 16 with the RHAl set of parameters. 

Fig. 6 The charge densities in 16 0, 40 Ca and 208 Pb. 

Fig. 7 The potentials of the nucleon and the anti-nucleon in 208 Pb. 

Fig . 8 The single-particle energies of the nucleon and the anti-nucleon as a function of 
density. The critical point for the NL1 is around 3.3 po while for the RHAl is around 
9.5 po- 
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Fig. 1 
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Fig. 4 
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Fig. 6 
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Fig. 8 




